﻿ Studiul transformatei Fourier discrete 1 Scopul lucrării Se studiază aplicaţiile ale TFD în analiza armonică a semnalelor şi în efectuarea unor operaţii liniare asupra acestora: convoluţia şi corelaţia Se utilizează în acest scop mediul de programare Matlab 2 Definiţii şi proprietăţi ale TFD Se consideră o secvenţă :xZC→ astfel încât [0xn= pentru nN Se ] 12 defineşte TFD a sa ca fiind secvenţa :XZC→ definită prin 22Nπ jk n− ⎧ N [] ⎪ ∑ xne N k N≤≤ 12, Xk= [] (1) = 1nN⎨ ⎪ 0, in rest ⎩ cu 1NN N= (2) 21−+ Uneori se prelungeşte prin periodicitate, de perioadă N dată de relaţia (2), semnalul [ xn ] si se obţine semnalul [ xn: ] ∞ [] xnxnmN=− (3) [ ] ∑ m=−∞ Notând: 2π j− New= N (4) se consideră semnalul: kn (5) [] [] ∑ NXkxnw= nN= semnifică faptul că însumarea se face pentru N unde în suma anterioară, notaţia nN= valori consecutive ale semnalului x[n] (o perioadă ) Se vede că semnalul [Xk este ] periodic, de perioadă N, ş i se poate arăta că: [ ≤≤ (6) ][] XkXk NkN= 12, Formulele de reconstrucţie ale semnalului iniţial din TFD a sa sunt : 2N 1 kn − [] (7) [] ∑ Nxn X kw= N= 1kN 1 kn− [] (8) [] ∑ Nxn X kw= kNN= În cazul utilizării ultimei relaţii rezultă: =−−− (9) [ ][][][] {} 12xnxn nN nNσσ s-a notat treapta unitate discretă) (cu [ nσ ] Se consideră, în continuare, secvenţele [xn cu transformatele corespunzătoare ] i principalele proprietăţi ale TFD, exprimate însă asupra [] iXk, i=1,2 Se prezintă semnalelor periodice asociate [ ][] ixn, respectiv iXk a) Liniaritatea [ + ∈ (10) ][][][] 12 1 2, , Cax n bx n aX k bX k a b+↔ b) Teorema întârzierii 0kn ↔∈ (11) [ ][] 00, Zxn n w X k n− c) Teorema deplasării 0kn− [ (12) −∈ ][] 00, Zwxn Xkk k↔ d) Teorema convoluţiei circulare ia de convoluţie circulară a secvenţelor periodice, de aceeaşi perioadă Notand cu operaţ N, al carei rezultat este secvenţă periodică, de aceeaşi perioadă, definită prin: − [] (13) [ ] [ ] ∑ 12xnxmxnm= mN= are loc: (14) [ ][][][] 12 12xnxn XkXk↔⋅ e) Teorema produsului 1 [] (15) [][] [] 12 1 2xnxn Xk Xk⋅↔ N f) Teorema lui Parseval 221 [] (16) xnXk= [] ∑∑ nN kNN== 3 Algoritmul TFR pentru calculul TFD Dezvoltarea prelucrării numerice a semnalelor a început odată cu elaborarea algoritmului de transformare Fourier rapidă (TFR, FFT conform iniţialelor din limba engleză) Acesta exploatează anumite proprietăţi ale exponenţialelor complexe Implementarea algoritmului TFR este eficientă când lungimea suportului secvenţei de transformat, respectiv a perioadei semnalului periodizat este o putere a lui 2 Când această cerinţă nu este îndeplinită, se prelungeşte intervalul din Z pe care se face transformarea până ce lungimea sa devine o putere a lui 2, completând semnalul cu eşantioane nule 4 Aplicaţiile TFD 4 1 Analiza armonică a semnalelor în timp continuu Această aplicaţie se bazează pe teorema eşantionării semnalelor de durată finită Se consideră un semnal în timp continuu, ( xt, astfel încât ) ( (17) supp 0,xtT⊂ ) [] {} şi semnalul rezultat prin repetarea prin periodicitate a acestuia (vezi fig 1) ∞ ()xtxtnT=− (18) ( ) T∑ n=−∞ Fig 1 Semnalul de analizat (a) şi repetarea sa prin periodicitate (b) Se notează cu () Xω transformata Fourier a semnalului (xt şi se presupune că (xt )) T admite dezvoltare în serie Fourier cu un număr finit de termeni: 2Nπ jkt T () (19) ∑ kTxt ae= kN=− Se eşantionează semnalul ( xt cu pasul: ) T T= E (20) 21N+ şi se consideră semnalul în timp discret: (21) [ () ] ExnxnT= având TFD X[k] În aceste condiţii se poate demonstra că, 12π ⎛⎞ [] (22) kaXk Xk== ⎜⎟ TT⎝⎠ adică elementele TFD sunt proporţionale cu eşantioanele transformării Fourier a semnalului () antioanele (21), xt În plus ()xt poate fi reconstituit cu exactitate din eş dacă pasul de eşantionare satisface (20) În realitate este puţin probabil ca această eşantionare să se facă astfel încât relaţiile de mai sus să fie îndeplinite, sau acest lucru este imposibil În acest caz, dezvoltarea lui () ine o infinitate de termeni: Txt conţ 2π jkt∞ T () (23) ∑ Tkxt ce= k=− ∞ şi în urma eşantionării se produce un fenomen de aliere a coeficienţilor: ∞∞ 12π ⎛⎞ ++ []Xk c X k N p== (24) () () ∑∑ ⎜⎟ (2 1)21kNp++ ⎝⎠ ppTT=− ∞ =−∞ Dacă valorile transformatei Frourier X(ω) sunt neglijabile pentru 2/NTωπ >, atunci valorile TX[k] estimează, în puncte k2π/T, k=-N, N, cu suficientă precizie eşantioanele acestei transformate În unele analizoare de spectru cu FFT se utilizează de obicei interpolarea (de obicei liniară) a eşantioanelor X[k] Dacă semnalul () lungă, se aplică tehnica „transformatei Fourier pe xt are o durată termen scurt” (fig 2) Prin această tehnică se analizează armonic, prin intermediul procedeului descris mai sus, porţuni din semnal, obţinute prin înmulţirea semnalului cu diverse ferestre Fereastra dreptunghiulară este cea mai simplă, dar, pentru evitarea apariţiei fenomenului lui Gibbs, se utilizează ferestre mai evoluate (Hamming, Gauss, etc), vezi fig 2 Fig 2 (a) Semnalul de analizat; (b) funcţia fereastra dreptunghiulară; (c) Semnalul analizat efectiv 4 2 Calculul convoluţiei Se consideră două secvenţe [ şi [, astfel încât: ]] 1xn2xn xnnNnN= [] 1120 pentru ş (25) xnnNnN=i <> [] 2340 pentru ş Pentru calculul convoluţiei celor două secvenţe se procedează în felul următor: se formează semnalele [xn şi [xn prin periodizarea semnalelor date cu o • ]] 12 perioadă −+ (26) 1NN N N N≥− 2143+ care să fie o putere a lui 2; în acest fel din convoluţia circulară a celor două secvenţe se poate recupera convoluţia obişnuită a semnalelor date • se calculează [] [] []XkXkXk= 12 se calculează TFD inversă, [xn a acestei secvenţe, reprezentând convoluţia • ] i [xn circulară a secvenţelor [ ]] 1xn ş2 • se calculează x[n] conform relaţiei: [[] [ ] [ 1]xn xn n N N n N Nσσ=⋅ −−−−−− (27) {} ] 13 2 4 4 3 Calculul corelaţiei Se consideră aceleaşi semnale ca şi la punctul anterior (ecuaţia (25)), dar considerate cu valori reale Funcţia lor de intercorelaţie se defineşte prin: ∞ =−∗ [] [ ] [ ] [ ] []Rnxmxmnxnxn=+ (28) 1212 1 2XX∑ m=− ∞ Pentru evaluarea acestei funcţii se procedează ca şi la punctul anterior, efectuându-se convoluţia între primul semnal refelectat în timp şi al doilea semnal 5 Desfăşurarea lucrării 5 1 Se consideră urmatoarele semnale: (a) () 1xt, din fig 3 a (impuls dreptunghiular) (b) () 2xt, din fig 3 b (impuls triunghiular) σ=−−⎡⎤ (impuls sinusoidal) (c) () ( ) ()() 30sinxtA t t tTωσ ⎣⎦ Fig 3 Semnalele ( )() 1xt si 2xt mân la alegerea utilizatorului programului „lfft” Se calculează Parametrii A, T şi ω 0 ră X, (ωX si (ωX şi se reprezintă grafic transformatele Fourier în timp continuu ( ))) 1ω23 Se calculează TFD ale semnalelor în timp discret ce rezultă prin eşantionarea semnalelor date începând de la t=-T/2 în 8, 16 si 32 de puncte pe durata T; se va utiliza în acest scop programul „lfft” Se compară ultimii doi membri care apar în relaţia (28) (vezi şi (24)) pentru fiecare eşantion în parte Se calculează TFD ca mai sus utilizând 256 de puncte de eşantionare Se compară figurile ce se obţin prin execuţia programului cu graficele transformatelor corespunzătoare şi se determină valoarea maximă a erorii 5 2 Se calculeaza convoluţiile: () , (xtxt∗, (xtxt∗ ())())() 11xtxt∗1222 ale semnalelor de la punctul 5 1 Se reprezintă grafic rezultatele Se eşantionează apoi în 512 puncte aceste semnale şi se calculează convoluţia lor, folosind programul “lfft” Se compară graficele calculate cu cele obţinute folosind mediul Matlab 6 Exerciţii în MATLAB În programul MATLAB pentru determinarea transformatei Fourier discrete se foloseşte funcţia fft Denumirea sa reprezintă prescurtarea de la Fast Fourier Transform (Transformata Fourier rapidă) şi indică faptul că este folosit pentru calcul un algoritm rapid Pentru obţinerea unor rezultate optime în cazul utilizării acestei funcţii este fie aleasă ca putere a lui 2 recomandat ca lungimea transformatei N să 6 1 Transformata Fourier în timp discret Să se calculeze transformata Fourier discretă a secvenţei: nπ ⎛⎞ [] 1sinxn= ⎜⎟ 5⎝⎠ Să se reprezinte grafic partea reală, partea imaginară, modulul şi faza transformatelor Fourier discrete calculate n=0:20; x2=sin(n*pi/5); X2=fft(x2,512); figure(1) subplot(2,1,1), plot(w,fftshift(abs(X2))),grid subplot(2,1,2), plot(w,fftshift(angle(X2))),grid figure(2) subplot(2,1,1), plot(w,fftshift(real(X2))),grid subplot(2,1,2), plot(w,fftshift(imag(X2))),grid Exerciţii Utilizând exemplul de mai sus, determinaţi şi reprezentaţi grafic partea reală, partea imaginară, modulul şi faza tranformatelor Fourier discrete pentru următoarele semnale: 1 [ pentru 010n≤≤ ] 1xnn= nnπ π jj−− 57 2 [] pentru 030n≤≤ 2xn e e=+ ,0 10nn≤≤ ⎧ 3 [] ⎨ 3xn= 20 , 11 20nn−≤ ≤ ⎩ nπ j− 52 ≤ 4 [] pentru 030n≤ 42xn e= 21nn++ 6 2 Relaţia de echivalenţă a energiilor 11NN−− 221 Pentru o secvenţă [ se demonstrează că: []xnXk= xn[] ] ∑∑ 00nkN== x=[(1:128) (128:-1:1)]; X=fft(x); a=sum(x *x) b=round(sum(abs(X) ^2)/256) Verificaţi cu ajutorul unui program MATLAB relaţia lui Parseval în cazul secvenţelor: = ,pentru 0,7nn⎧ ⎪ 1 [] ⎨ 1xn= 0, în rest⎪ ⎩ 4, pentru n=0⎧ ⎪ 2, pentru n=1 ⎪ ⎪ 2 [] ⎨ 21, pentru n=2xn= ⎪ -1, pentru n=3 ⎪ 0, în rest⎪ ⎩ 